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ABSTRACT 

ANISOTROPIC SOLAR COSMIC RAY PROPAGATION 
IN AN INHOMOGENEOUS MEDIUM 

A phonomenological , 1-dimensional model of solar cosmic ray 
propagation is presented which allows one to compute both the 
anisotropy and the intensity as a function of time during the early 
phases of certain events. It applies for arbitrarily large 
anisotropies, ranging from 0 to 1, £ind is valid for even the earliest 
times. Thus, it is not subject to the restrictions of earlier 
models based on the diffusion equation of telegraph equation. 

The model assumes that cosmic rays are released instantaneously 
and effectively interact with isolated scattering centers which are 
distributed uniformly throughout the solar wind. Each scattering 
center is assigned a probability P(x^) that a cosmic ray will be 
reflected by it; this function ranges from 0 to 1 and can be any 
function of position. The mathematical formulation of the model 
consists of a set of difference equations involving P(x^) and the 
probabilities f^ and f~ of finding a particle moving away from or 
toward the source, respectively, at a given time and position. 

The model is used to examine the general effects of a diffusing 

region close to the sun (the solar envelope) on solar cosmic ray 

propagation. The model, applied to the >1.1 bv protons of the May 4, 

1960, event, gives a good fit to the intensity- time curve for the 

"sunward" flux and correctly gives a negligible "anti-sun" flux through 

out the event. The characteristic dimension and diffusion coefficient 

21 2 , 

of the diffusing region are found to be 27 R© and D = 10 cm /sec, 

max 


respectively. It is suggested that the solar envelope may be identified 
with an extended solar wind heating region and that the cosmic-ray 
scattering may be due to the hydromagnetic waves which Barnes has 
postulated as the energy source. 
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ANISOTROPIC COSMIC RAY PROPAGATION IN AN INHOMOGENEOUS MEDIUM 

lo The Solar Envelope 

L. F. Burlaga 

I. Introduction 

The existence of a diffusing region close to the sun has been inferred 
from recent observations of low energy, solar cosmic ray anisotropies 
(McCracken et alo, 1967 , Lin et alo (1968)) and from an early study of the 
anisotropy and dispersion effect in a high energy (>1 Bev) solar cosmic 
ray event (Lust and Simpson, 1957). Following Lust and Simpson, we 
shall call this region the "solar envelope". The above authors, also infer 
the existence of a second diffusing region near the earth and beyond, 
which is the cause of the trend toward isotropy with increasing time. 
Following Parker, who was the first to study this region quantitatively, 
we shall refer to it as the "solar shell." The large anisotropies observed 
during the rise of certain events imply that relatively little diffusion 
occurs between the solar envelope and the solar shell. 

There are many diffusion madels which attempt to examine various 
aspects of the propagation, but as yet there is no model which properly 
accounts for the highly inhomogeneous structure of the solar wind as 
described above, and which describes both the solar cosmic ray intensity 
and anisotropy as a function of time. The aim of this paper is to present 
a simple, phenomenological, 1 -dimensional model, valid even in the limit 
of large anisotropies, which enables one to use the measured anisotropy 
and intensity profiles of solar cosmic rays to study the inhomogeneous 


structure of the medium. 


V n 
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The model is presented in Section 2<« In Section 3 the model is 
applied to examine the general effects of the solar envelope on cosmic 
ray propagation, and it is shown how the model can be used to analyze 
certain types of cosmic ray events o The unusual event on May 4, 1960, 
is analyzed in Section 4, where it is shown that the model gives good 
fits to the intensity- time profiles for both the forward and baclcward 
fluxes, and that the size of the solar envelope is '^27 Rq. A possible 
relation between the solar envelope and the solar wind heating region is 
also discussed in Section 4« 

For most events, the solar shell must also be considered, because 
this is the cause of the approach to isotropy.. This will be considered 
in a subsequent publication, where analytical properties of the model 
and its relation to models based on the Boltzman equation will also 


be discussed 
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II. The Model 

We consider propagation in a 1 -dimensional, semi-infinite medium 

with a source at xs=0. The medium is represented by an infinite 

number Of point scattering centers which are equally spaced and extend 

from the source to infinity. The scattering properties of the medium 

are described by two parameters: K, the constant separation between the 

scattering centers, and the "effectiveness" of the scattering centers, 

P , which is the probability that a particle will be reflected when it 
L 

encounters a soatterlns center at The properties of the cosmic 

rays are represented by three parameters: a) their speed v, 

which is assumed to be constant, b) f , the probability that a 

L, I 

particle is moving outward (i.e. away from the source) after a collision 

at L-1, T-1 and is just approaching point L at time T and c) 

the probability that a particle is moving inward (i.e., toward the 


source) after a collision at L+1, T-1 and is just approaching point 
L at time T. Thus, f is proportional to the flux of particles 

Xj iL 

moving away from the source (the forward flux) and £ is proportional 

L,T 

to the flux of particles moving toward the source (the backward flux). 


One can write the following recurrence equations for f"^ and f“: 

r* = ? *0-rA-f*T 

T + e^T+l ^ hT ^ J L^T 




= (i-r,)i^T Ct 


The first of these equations states that the probability that a 
particle is moving outward toward L+i at time T+1 is the sum of the 
probability t!<at a particle moves inward toward L at T and is 


m m 


reflected at that point, plus the probability that, a particles moves 
outward toward L at T and is not reflected at that point# The 
second equation has a similar meaning. These equations form the basis 
for the model. Given a source function, and the boundary conditions 
one can use (1) to calculate f ^ L T 

For the source function, we choose 


r s 1; 
0,0 


f « 0, 
0,0 


( 2 ) 


corresponding to an instantaneous point source. For the boundary 
condition, we choose 


- f 

^ 1,T o,T>l’ 

which implies a reflecting boundary at x =0 

L 


( 3 ) 


• Application ty fchti golai ; 

Idually, ono v;ould like to uae tiiu model to deduce F(‘x) . Of 
course, this cannot be done. However, one can proceed inductively 
and use the model to test certain hypotheses concurning F(x) . This 
section follows the second approach and eKamine?? the effects of a 
diffusing region close to the sun (the solar envelope) on solar cosmic 
ray propagation. 

Let us assume that 

I 2 

F(L) = .5e (4) 

The choice of a gaussian is somewhat arbitrary, but the results of 
interest do not depend sensitively on the exact form of P(L) . The 
coefficient .5 was chosen because in the limit the model 

reduces to the 1 -dimensional random walk in a homogeneous medium, if 
P » .5. The essential effects of the solar envelope on cosmic rays 
are determined by the number of scattering centers in the envelope, 

which is measured by L . 

D 

To illustrate the properties of the model and the effects of the 
envelope, we have calculated the sunward and antisunward "fluxes", 
r and f , for various values of Lj^, assuming that there are L^ 
scattering centers betv/een the observer and the source, The results 
fall into three classes, corresponding to L «L , and L J5L .. 

when L_<<1, , f“ is found to be negligible, as one would expect, 

D o 

since there is a negligible probability of backscattering beyond the 
observer's position. The characteristics of f*^ are illustrated in Figure la. 
When Lp<2, the maximum intensity occurs at the onset, Tq, and ft*T) decreases 


•• 6 - 

exponential3y at all times. When Lj)= 6 , the envelope begins to act 
like a diffusing region; the particle intensity rises to a maximum 
at Tm, shortly after onset and then decays exponentially as result 
of the loss of particles from the 'diffusing' region. As Lj^ increases 
to 10 , the maximum shifts to later times and the decay, though still 
exponential, is less rapid as a result of the greater effectiveness 
of the diffusing region in delaying the escape of particles. With 
the observer at L^=50, the values Lj^ = 2,6, and 10 correspond to a 
solar envelope with dimensions 8.4 R©» 25.2 R©, and 42 R 0 , respectively. 

When the scattering centers beyond Lq effective, 

and an appreciable anti sun flux, f is observed,. This is illustrated 

in Figure lb for and Ljj=40. The rise to maximum and exponential 

+ - 

decay are seen in both f and f , and the trend toward increasing 
Tmax with increasing Lq continues o The anisot-.'opy approaches 

a limiting value which decreases as increases. This is 28% for Lj^= 
40 and 68 % for Lj^=30o 

When Lj^»Lq 5 the propagation tend,s to be like 1-dimensional 
diffusion, as one would expect, so this' cas® will not be discussed 
further. 

Having illustrated the general effects of a solar envelope of 
various dimensions on cosmic ray propagation, we now show how certain 
cosmic ray observations can be used to infer the size R and diffusion 
coefficient D = Xv/2, of the solar envelope, on the assumption that 
P(L) has the general form of (4) . 
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Since 


and 



2 


D(— ) 
'sec' 


x210 (in units of R©) 

_ Io5xl0^\ (cm/sec) 

2 L 


(5) 

( 6 ) 


o 

the problem can be stated, "Given f (t) and f (t), and assuming P(x) 
is given by (4), find and Lp.’’ 

From numerical calculations with ranging from 50 to 500 and 
Ljj ranging from 5 to 50, the following empirical relation was found 


L 

D 


?=100 


(T^ / Iq-1)^ 

<’^D ! 


< L^/4, 


(7) 


The error in this relation is ies^ than 10%. When f’*'(t) 
has the form shown in Figure 1, it can be characterized by three 
parameters: 1) the onset time, t^ * x^/v, where x^ is the position of 

the observer, 2) the time of maximum, t^jj, and 3) the decay time tj^. 
Since t = Tt , where t = \/v, we have T_/T^ = t /t and T_/ T = ^ /t . 
Thus, can be determined Very easily from the observed using (7). 

Similarly, it was found that there is a simple relation between 
Tp and Lj^, which is given in Figure 2» This is independent of when 
Ljj<L^, indicating that the observer sees the same decay profile at 
any position beyond the "diffusing" region and that the decay constant 
is determined only by the number of effective (P'^o5) scattering centers 
in the envelopeo Thus, given Lj) from (7), one obtains Tj) from Figure 2. 
Lq can then be determined from the relation 


L 


o 





( 8 ) 


Thus and are determined by (7) and (8) if t^, and tjj are 

known, qoOcd., R and D are then given by (5) and (6). 
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IV . The May 4, 1960 Even t 

We ahall ii show how the model can be applied to give the 
characteristics of the envelope and the intensity -time profiles of 
the forward (looking toward the source) and backward (looking away 
from the source) fluxes for the >1 bv protons in the May 4, 1960, 
event. 

The classical May 4 event presents a special challenge to propagation 
theories, for it was of unusually short duration, the anisotropy was always 
essentially 100%, and it defies any attempt to describe it by the 
common diffusion models. On 'the other hand, it was clearly associated 
with an optical flare and a cm radio burst at 1025 UT, so that the 
injection time can be determined (assuming instantaneous injection); 
and good measurements of the intensity- time profiles for several 
viewing directions are available, so that f'’’(t) and f"(t) are well- 
known. Moreover, the flare occurred at 80°W, which is only ~10° 
from the theoretical position of the base of the interplanetary 
magnetic field line that passes through the earth (Burlaga, 1967), 
so the use of a 1-dimensional model is justified if perpendicular 
diffusion can be neglected near the sun as suggested by Axford (1965), 
Burlaga (1967) and others » 

The data which we use for this event are taken from McCracken 
(1962), Figure 2.2<. For the forward flux, we use the data from the 
Fort Churchill neutron monitor which was viewing at an angle 6''40*^ 
from the streaming direction. These data are shown by the dots in 
Figure 3. The backward flux is inferred to be zero for all t, because 
Mirny, which viewed at 6'^120° saw no enhancement at any time. Since 


. 9 ~ 

£ was aXvmys near the background level, and since was characterized 
by a rise to maximum followed by an exponential decay, it is likely 
that most of the scattering occurred near the sun, so that the model 
of the preceding section should be applicable for this particular 
event. 

From Figure 3 wo find that t,^ = 13.3 min, t =22.5 min. and 
« o ' m 

tjj = 20,1 m in. Thus, (7 ) gives = 20.9, Figure 2 gives T^ = 241, 

and ( 8 ) gives = 160, The collision probability P(x) (equation 

( 4 )) is now determined, and the model can be used to calculate 

f and f . The results for f (t) are shewn by the curve in Figure 3, 

which is normalized to the maximum intensity. It can be seen that 

the theory gives an excellent fit to the observations. The corresponding 

results for f"(t) (which, recall, is on an equal footing with f^(t) 

-24 

in our model) show that the maximum is <10 ^ in the units of Figure 3. 

Thus the predicted backward flux is essentially zero, in agreement 
with observations. We conclude that the model provides a satisfactory 
descript i.on of the observed forward and backward fluxes of the 
>1.1 bv protons in the May 4, 1960 event. 

Now let us examine what the model tells about the solar envelope 
on May 4, 1960. Since = 160, implying 160 scattering centers 
between the sun and the earth, the separation of the scattering centers 
was X = .00625 ATJ. This is approximately equal to the mean free path 
near the sun. Using ( 6 ) with v=lo9xl0 cm/sec., we find that 

the maximum value of the diffusion coefficient in the solar envelope 
was ]>^xl0^®cm^/sec. Similarly, using = 20,9 we find from 
( 5 ) that the size of the diffusing region was R=27 Rq. 
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Note that the decay time tjj « 20 mino , is approximately the time required 

2 

to diffuse out of the shell, R /(4D) = 17 min. 

Our diffusion coefficient may be compared with the estimate of 

20 2 

Shishov (1966), D^^SxlO cm /sec, based on the assumption that Rr-IOR 

0 

on May 4, I960, His estimate is inferred from an analytical model for 
propagation in a spherically symmetrical, homogeneous medium from a 
source of the form exp(-t/16 min). While our model confirms Shishov' s 
intuition, it differs fundamentally from his approach and provides more 
precise results. Our model might also be compared with those of 
Axford (1965) and Fisk and Axford (1969). They presented equations for 
f^ and f , but these are for a homogeneous medium. The solutions cannot 
describe the May 4, event, We also note that Reid (1964) and Axford 
(1965) presented models, postulating a very thin (R« 1 R 0 ) envelope to 
give diffusion around the sun. 

Finally, let us consider the origin of the diffusing region, 

Burlaga and Ogilvie (1969) presented evidence in support of the hypothesis 
that there is an extended heat source near the sun which accelerates 
and heats the solar wind, Hartle and Barnes (1969) have shown the 
relation between bulk speed and temperature of the solar wind, found 
by Burlaga and Ogilvie, can be obtained by assuming that protons are 
heated out to ^^20 Rq at quiet times. Note that the size of the heating 
region is approximately the size of the solar envelope inferred from 
our model. Thus, we offer the hypothesis that the heating region may 
be identified with the solar envelope. This is supported by the 
suggestion of Barnes (1969) that the heating is produced by damping 
of hydromagnetic waves, for these waves cou’ ’ act as cosmic ray 
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scattering centers. This hypothesis raises the interesting possibility 
of using cosmic ray observations and our model to study the solar 
wind heating region. Since Hartle and Barnes find that a larger 
heating region is required to give higher solar wind temperatures, 
we would expect that the size of the envelope inferred from cosmic 
ray observations would be positively correlated with the solar wind 


temperature 0 
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FIGURE CAFTIONS 

Variation of f*^ and f" with T for a particle released 
at T=0 into a medium with P(L) » The 

observer is at Lq* 50. For < 10 (panel (a)), f” is 
negligible. As approaches L^, f" approaches f'*' 
(panel (b)). 

This shows how Tp is related to L The relation is 
Independent of L^, when Lj^ < h^/4. 

The points are observations of >1.1 Bv protons seen 

4* 

at Fort Churchill on May 4, 1960. The curve shows f 
obtained from the model used to generate the curves in 
Figure Ij no free parameters were used to get this fit. 
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